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1 Introduction 



Fractional Brownian motion is a one parameter family of zero mean Gaussian processes 
indexed by H £ (0, 1) which was mathematically introduced by B.B Mandelbrot and 
J.W. Van Ness in [5], (cf. [I] as well). It generalizes Brownian motion by setting 
EBf = t 2H , and can be used to model various phenomena, in finance as well as in other 
fields. This is primarily due to the fact that its self-similarity depends on the parameter 
H (named after H. Hurst, [3]), which allows for phenomena exhibiting different kinds 
of self-similarity to be modeled by one of the fractional Brownian motions. 

Since fractional Brownian motion is not a semimartingale (unless H=^), the ordi- 
nary stochastic calculus for semimartingales (such as the Ito integral) does not apply. 
Instead, there are several approaches for defining a stochastic integral with respect 
to fractional Brownian motion. The divergence integral is one possible approach, us- 
ing the Malliavin divergence operator as the basis for integration, a survey of which 
can be found in D. Nualart's book [5J. A different approach was developed by Zahlc 
in [8] which involves a pathwise definition of the stochastic integral. This requires a 
generalization of the Young-Stieltjes integral, introduced in the same paper. 

Given a suitable process u, the divergence integral t h-> J Q u T dB T yields a new pro- 
cess X t . While the general theory of Malliavin calculus provides an abstract formula 
for the covariance function of X t (cf. [6j), in many concrete cases it is not straightfor- 
ward to express it explicitly. This paper provides an (integral) such expression when 
u T = F(B T ), for an appropriate function F. More generally, if X t = J Q F (B T ) dB T 
and Y s = G (B<j) dB ai it will be shown in Proposition 13.11 that 

EX t Y s = f f (f (r, a) + /i (r, a)) dr da (1.1) 
Jo Jo 

where /o (resp. fi) is an explicit kernel involving F and G (resp. F' and G'). This 
formula will be later extended in Theorem 13.61 to non-smooth F, G. 

In [2] Y. Hu and D. Nualart showed that if H ^ 1 /2 the process X t = f Q sgn(B s ) dBs 
is not a fBm (unlike the Brownian motion case). In the absence of a formula for the 
covariance of this process, a detailed analysis of its chaos expansion was necessary to 
reach that conclusion, and this was one of our motivations to obtain such a formula 
which, moreover, would indeed have to accommodate non-smooth functions F such as 
sgn(x). 

Section [2] includes some preliminaries and two auxiliary results fLemma l2.2[ whose 
proofs is deferred to the Appendix, and Proposition I2.3|) . The main covariance for- 
mula (jl.lj) is presented, with explicit expressions for the kernels /o and /i, in Subsec- 
tion [3T] for smooth functions F and G (Proposition 13. 1|) and in Subsection 13.21 where 
the regularity assumptions on F and G are significantly relaxed (Theorem 13. 6[) . 

This work is based on the first author's Master thesis. 
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2 Preliminaries: Fractional Brownian Motion 



The following introduction to fractional Brownian motion and its analysis is taken 
mostly from Chapter 5 in [6]. 

Definition 2.1. For H £ (0, 1), a process {B t , t £ [0,T]} is a fractional Brownian 
motion with Hurst parameter H if it is a centered Gaussian process with covariance 
function 

R H (t,s) = ±(t 2H + S 2H -\t- S \ 2H ). (2.1) 

The parameter H will not be explicitly indicated in the notation Bt. It follows that 

E\B t -B s \ 2 = \t-s\ 2H (2.2) 

and, by Kohnogorov's continuity criterion, we may assume that {B t | t £ [0, T]} has 
a-H61der continuous trajectories for any a < H. When H = 1 /2, B t is a standard 
Brownian motion. 

From this point on, it will be assumed that H > i . The discussion takes a different 
turn when H < i (more about that in Section!?]). It will be convenient to write 



R}i{t,s)=a.H I I |t — <r| 2ir 2 drdcr, 
Jo Jo 



(2.3) 



with 

a H = H(2H-l), (2.4) 

as can be easily verified. Denote by £ the space of step functions (that is, the space 
spanned by all indicators of subintervals of [0,T]), and on it define 

(f,9)=a H [[ f{T)g{o)\T~o-\ 2H - 2 &T&o-, f,g££. (2.5) 
J J [o T ] 2 

In particular, (l[o,^, l[o,*]) = Rif(t,s). The completion H of £ with respect to (•,•) 
plays a fundamental role. It contains distributions as well as proper functions, but the 
function space 

\H\ = j/: [0,T]->M | ||/||f K| :=a H J ^ \f (t)\ \f (s)\ \t - s\ 2H ~ 2 dsdt < ooj 

(2.6) 

is a (strict; cf. )7]) dense subspace of H, while the inner product of H on elements of 
\H\ is still given by 

The linear isometry l\o,t] 1 — ► <^l[o.t] := B t £L 2 (P) extends to H. The Malliavin 
gradient D is then defined as an unbounded operator from L 2 (P) to L 2 {P;'H) by 
setting DSh = h and extending it, by the chain rule 

n 

f (Shi, . . . , Sh n ) i — > ^2 ®if (^li • • • > Sh n) h i 

i=l 
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and closure, to its domain B 1 ' 2 (the same can be done for any p > 1 instead of p = 2). 
The Malliavin gradient of F G I? (P; K), for a separable Hilbert space K, is defined 
similarly as an element of L 2 (P; Ti (g> if), its domain denoted D 1 ' 2 ^ ) accordingly. We 
now state a lemma which relates the differentiability of a K- valued variable to that 
of its linear marginals. This auxiliary result is quite standard, but for the sake of 
completeness its proof will be included in the Appendix. 

Lemma 2.2. Let u G L 2 (P, K). Then u G B 1 ' 2 (K) if and only if for each k G K we 
have (u,k) G D 1,2 and there exists an orthonormal basis {e n }^Li in K such that 

oo 

J2E\\D(u,e n )f <oo. 

71=1 

In that case, Du viewed as a Hilbert- Schmidt operator from K to Ti is given by 
Du(k) =D(u,k), k E K. 

The dual operator Su acting on suitable elements of L 2 (P; Ti) is commonly referred 
to as the stochastic integral J u t dB t , namely 

E (DF, u) = EFSu VF G D 1,2 . 

The domain of <5, Dom (6), contains D 1,2 (T-L) (cf. [HJ Proposition 1.3.1]), and moreover 

ESuSv = E(u,v) + E trace (Duo Dv) , u, v G B 1,2 (U) . (2.7) 

If u = h G % the notation 5u is consistent with the previous usage of 8h, and the 
term "integral" just reflects the fact that, for H — 1 /2, 8u coincides with the Skorohod 
integral of Ut with respect to Brownian motion (although when H ^ 1/2 there is no 
integral involved any more). 

Let u be in Dom(<5). In order to make sure that wl[o,t] too belongs to <5's domain 
for any t (recall that u(uj) is not necessarily a proper path), further requirements on 
u will be imposed following 6, p. 288]. Namely we shall assume u to be |?^| -valued 
and Du to be representable by an appropriate kernel in the following sense. Define 
\%\ ® \H\ to be the space of measurable functions K (s, t) on [0, T] 2 such that 

ll*H|tt|®|«| :=a H [[[[ \K(s,t)\\K(u,v)\\s-u\ 2H - 2 \t-v\ 2H ~ 2 dsdtdudv <oo. 

This forms a Banach space which, on the other hand, and equipped with the inner 
product 

(K 1 ,K2} H<sn = a H /// / Kx(s,t) K 2 (u,v)\s - u\ \t — v\ dsdtdudv, 

[0,T] 4 

can be isometrically embedded as a subspace of T-L® T-L (of Hilbert Schmidt operators 
which possess a kernel for their action on \H\). Thus, the subspace of D 1 ' 2 (TV) we shall 
assume u to belong to is 

B 1 ' 2 (\H\) = {ueB 1 ' 2 (H) I u(co)e\H\ a.s. and Du(w) G \H\®\H\ a.s.} . (2.8) 
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It turns out that whenever u belongs to D 1 ' 2 so does ul[o,t] for every t£ [0,T], 

and in particular 

f u s dB s = 5(ul [0 . t] ) (2.9) 
Jo 

is a well defined process. This divergence integral for fractional Brownian motion has 
corresponding ltd and Tanaka formulae. 

The last auxiliary result of this introductory section, possibly interesting in its own 
right, concerns a class of Hilbcrt-Schmidt operators on W, to which Du belongs - as 
suggested above - when weB 1,2 {\H\) and certain boundedness conditions are imposed 
on its gradient. 

Proposition 2.3. Assume that the HUbert- Schmidt operators K\ and K 2 on H are 

characterized by bounded kernels k\ and k2, in the sense that 

K l h = a H dt h(t)\t- s\ 2H ~ 2 ki(s,-) ds, he\H\, €=1,2, 
Jo Jo 



where au is defined in (12. 4p . Then the trace-class operator K — K\K 2 is characterized 
(in the same sense) by the kernel 

k(s,t)=a H dr k 1 (s,T)k 2 (a,t)\T - <t\ 2H ~ 2 dcr. (2.10) 
Jo Jo 

and its trace is given by 

trace (#)= a H / ds k {s, t) \s - t\ 2H ~ 2 dt. (2.11) 
Jo Jo 

Proof. The first claim involving (|2.10|) is a straightforward application of Fubini's 
theorem. As for the trace, recall that 

where the inner product is the Hilbert Schmidt inner product. Moreover, K\ is given 
by k\ (s, t) — ki (t, s). It therefore suffices to prove that for any Ki, K 2 satisfying the 
assumptions of the Proposition, 

(Ki,K 3 ) =a„ jk x {s u t x ) k 2 (s 2 ,t 2 ) \ Sl - s 2 \ 2R - 2 \h - t2\ 2H - 2 d Sl ds 2 dhdt2. (2.12) 

[0,T]-» 

When Ki = h\ ® h\ where hj € |%| and i,j — 1, 2, (|2.12ll is a simple calculation. This 
clearly carries over to sums of such operators. For general k\ and k 2 , let {k™}^ =l be a 
sequence of kernels of the type above which converges almost everywhere to ki and is 
uniformly bounded (this sequence exists since fc, itself is bounded) . Then each kernel 
defines a Hilbert Schmidt operator if™, and (|2.12j) is known for K™. It remains to 
prove that the sequence {K^}^ =1 converges in the Hilbert Schmidt norm to Ki, and 
that the right hand side of (|2.12p converges accordingly. 
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To see that K™ n -Z^>? note that {K™} 00 . is a Cauchy sequence, as implied by 
dominated convergence with the already established formula in (|2.12|) for the Hilbcrt 
Schmidt norm \\K™ - Kf\\. For tmyhe\H\, 

K?h (t) = a H [[ h (s) fc™ (r, t) \s - t\ 2H ~ 2 ds dr 



'[0,T] 2 

which converges in H to Kih by dominated convergence. Thus ™^°p Thus 

lim (K?,K?) = (K 1 ,K 2 ). 

n— >oo 

To see that 

/2^^ 2 2 
fc™ (si,*i) k~2 \S2,t2) \si - S2I |*i — *2] dsids 2 d*idt 2 
[0,T] 4 

fci (si,ti) fc 2 {s2,t 2 ) I si - s 2 | |ti — *2| dsids 2 d*idt 2 , 

[0,T] 4 

note that the integrand on the left hand side is bounded by some constant multiple 
of the integrable function |si — s 2 | |*i — ^2! • An additional application of 
dominated convergence concludes the proof. □ 



3 The Cross- Covariance Formula 

We recall from Section [2] that a process u in B 1,2 belongs to the domain of 8 and that 
for any u,v£V) 1,2 , formula (|2.7[) holds, namely 



E I u T dB T / w r d£> T = i? (u, u) + E trace (Dm o Z)?;) 
Jo Jo 



(3.1) 



where Du is viewed as a Hilbert-Schmidt operator-valued random variable. Our aim 
is to find a concrete expression for the right hand side of (|3.1[) when the integrands 
are respectively of the form 

u = F{B.)l m {-) and v = G{B.) l [o , s] (0, (3.2) 

for s,i>0, initially for sufficiently smooth functions F and G. 

3.1 Regular Coefficients 

Recall the constant ctH=H(2H—l) which was already defined in (|2.4p . 

Proposition 3.1. Let F and G be twice continuously differentiable functions on R 
with bounded derivatives of first and second order, and let u = F (B.) and v = G (B.). 
Then u,v £ B 1,2 (|%|) (defined in (I2.8[) ) and the following formula holds: 

drder 



E [F{B T )dB T [G{B a )dB a =a H [ f f|r - a\ 2H 2 M (r, a)+ 7 (r, a) P (r, a) 
Jo Jo Jo Jo L 



(3.3) 



G 



where M (r, <j) =EF (B T ) G (B a ) , P (r, cr) =EF' (B T ) G' (B a ) and 
if 



7 (t,ct) 



2il- 1 



/ , , 1 i \2H—1 

(t V cr) — \t — cr | 



(TAaf^ + lr-fTl 2 ^ 1 



(3.4) 



(The first term in the right hand side of (|3.3[) corresponds to the isometric component, 
familiar from the Brownian setting, and the second term could be seen as a correction 
when H> 1 /2.) 

The proof of Proposition 13.11 will consist in identifying Du and Dv, on the right- 
hand side of formula (|3.1I) . as kernel operators (Lemma I3.2[) . and then the trace of 
Du o Dv with the help of Proposition 12. 31 



Lemma 3.2. Assume that F 6 C 2 (R) /ias bounded first and second derivatives, and 
let u t — F (B t ). Then u £ D 1 * 2 (\H\) and Du, viewed as a Hilbert- Schmidt operator, 
acts on [H\ according to 



(Du)h(-) =a H I ds I h (s) \s - t\ 2H ~ 2 l [0<t] (•) F' (B t ) dt 



(3.5) 



Proof of Lemma [KR Note that, almost surely, u is bounded and thus belongs to \%\. 
Denoting tk — = the smooth approximations u n = Y^k=i F l[t fc _i,t fc ] 
converge to u in L 2 (P, %). This can be easily seen by writing 



u t - u 



Next, u n —u n+1 = J2ti [J , (S.«- a )-- F ( B ««-i)] ifc-,*], where Sj =^ +1 = 
so that 



T 



2 if 2" 



IV - ^«+i 2 < 2" — T J2W F ' ( B ^) Mo 



The terms in the sum can be estimated by 



'-jj-a] ~ F ' (^2,-1) l[0,s 2i -i]| 



l-P 1 ' (Aisi-a) 1 [0, S2 ,-2] ~ F ' ( B *ai-i) ilO^ai-i] || 2 

< 2 f [F' (S S2i _ 2 ) - F' (B^J] 2 | S2l _ 2 | 2ff + F' (B^) 2 {^) 2H 



whence 



E\\Du n -Du n+1 \\ < 



1-iH ( ll^//|| 2 j,2H 



\F'\\t)T 



^4H 



(2 4H - 



This implies that {Du n }^ =1 is a Cauchy sequence in L 2 (P,ri €5 H) (H> 1 /2), so that 

u e p 1 ' 2 {%). 
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To obtain the expression (|3.5p . let ft £ \H\. Then 

2" 

(Du n ) h(t) = ^ (ft, l [tfc _ lltfc] > F' (B tk _J lp^^j (t) 

fc=l 

where 



#(r) = / ft(a)|a-r| 2J ^ 2 da. 
Jo 



It remains to show that the right hand side converges to the right hand side of (|3.5I) as 
n — > oo. But this is a standard application of the dominated convergence theorem. □ 

Remark 3.3. In fact, for each t G [0,T], a straightforward modification of the proof 
shows that D (ulfo^i) is given on elements of \H\ by 

D (ul m ) ft(-) = a H f da [ h (a) \a - r| 2 ^ 2 l [0 , r] (■) F' (B T ) dr. 
Jo Jo 

so that its kernel, in the sense described in Provosition \2.3l is 

h (t, a) = F' (B a ) l [Tjt] (a) l [0>t] (r) . (3.6) 

Proof of Proposition ^. 1\ Fix s,t € [0, T]. By Lemma l3~2l and its subsequent Remark, 
D (itlr ot i) is represented by the kernel (|3.6[) and -D (ulr 0)S iJ is similarly represented by 

k 2 (r, a) = G' (B a ) l [r , t] (a) l [0>t] (r) . (3.7) 

We use these kernels for the derivatives in (|3.1[) , and Proposition 12.31 to evaluate the 
trace: 
ft 



E Qf F (-B r ) dB r J G (Ba-) dB^j = 



E la B II F(B T )l m ( T )G(Ba)l l0 , s] (a)\r-a\ 2H 2 da At 



[0,T]- 



+E\a H I oh F' (B T2 ) l [Tlit] (r 2 ) l [0)t] (n) G' (S CT2 ) 1 [cti>s] (tr a ) 

T[o, s ] |r 2 - cri| 2ff_2 dT 2 d(Ti^ |t x - <J 2 | 2fl ~ 2 der 2 dn 

= a ff / dr / F(F( J B T )G(S CT ))|r-a| 2ff - 2 da + 
Jo Jo 

+a% f dn f da, f dr 2 fE(F' (B T2 ) G' (B a2 )) |r 2 - a^ 2 \ n - a 2 \ 2 "- 2 da 2 . 

Jo Jo Jt x Jax 

and the expression in (|3 . 3[) now follows by changing the order of integration (by where 
the integrals according to t% and a\ are carried out first, resulting in 7(t 2 ,<7 2 )). □ 
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3.2 Non-regular Coefficients 

Proposition 13.11 is unsatisfactory in that the smoothness assumptions on the func- 
tions F and G are too restrictive. In this section it will be extended to accommodate 
(equivalence classes of) real functions whose distributional derivatives can be decom- 
posed into an absolutely continuous component with bounded density with respect to 
Lebesgue measure, and a finite discrete component (whose support will be assumed 
for simplicity to have no accumulation point). Namely, denote 

ac = {F : E -> R is absolutely continuous, F' € L°°(M)} , 
Qd = { F = y^x n l( a a +1 ), y^ \x n+1 ~x„ \<oo, lim a n = ±00 > 

k n£Z nGZ ) 

and 

Q = | equivalence classes of elements in C? ac + under a.e. equality | 

Remark 3.4. The distributional derivative of any F € Q is a measure F'(dx). If, 
however, F has no discrete component, and by a slight abuse of notation, we shall also 
write F'{dx) = F'{x)dx. 

Lemma 3.5. Let F E Q. Then: 

1. L^RjletynCoCR) Ci^R.I-F'l) . 

2. F= lim F n pointwise, where each F n can be decomposed as the sum of two C°°- 

functions with bounded first and second derivatives, F n = F^ + F%, such that 
the sequences F^ and F% are uniformly bounded respectively in L°°(R) and in 
L (R) (thus the F n 's satisfy 1) ). In particular, there exists an M > such that 

\F (x)\ + \F n (x)\ < M (I + \x\) VieR.VrieN. (3.8) 

Moreover 

lim / hF^dx = I hF'(dx), WiG-L^R) n C (R) . (3.9) 

Proof. Assume without loss of generality that F is an increasing function and that 
either F e Q &c or F e Gd- 

1. If F e § ac , denote /= 4gg) . Then J R \h\ F> {dx) =f R \h\ fdx < WfW^ 11% < 
00. The case F £ G d follows similarly. 

2. Let {^n}^i be a symmetric approximation to the identity (that is, ipi <G (R) 
is an even non- negative function such that J R (p\ dx = 1 and ip n (x) — nipi (nx)), 
and define F n — F * (p n . Then F is the pointwise limit of the F n s. Moreover, if 
F e Oac then 

K\ = \F>*ip n \<\\F'\\ 00 , 
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whereas if F £ Gd then F' = Y^nLi ^ x i&ai {8 a denotes the Dirac measure at a) 
and 

„ oo oo 

ll-^nlli - / Axjfn (x - Oj) da; = Ctj < oo, 

JR i=l i=l 

Finally, if F £ Q ac then F' is the almost everywhere limit of F' n = F' * ip n 
and, since F^ — F' is bounded, the limit follows from the dominated convergence 
theorem. If, on the other hand, F £ Qd, then 

• 00. oo 

/ hF' n dx — y~] Axi / h (x) ip n (x - a») dx = V" Ax t h * cp n (a t ) . 
Jr i=1 i=1 

The continuity of h implies that lim„_ ! . oc h * tp n = h at every point. Since 
\Axih * tpn (ai)\ < Axi WhW^, 

p oo oo 

lim / hF'dx = lim > Axj/i * ip n (ai) — > aih(xi). □ 

n— »oo /to, n— >oo ^ — ' — ' 
JK i=l i=l 

Theorem 3.6. Let B t ,t £ [0,T], be a fractional Brownian motion with H > 1 /2 and 
F,G £ Q (or more generally satisfy the conclusions of Lemma \3.5\) . Then, almost 
surely, the processes u := F (B.) l[o,t] and v s := G (B.) l[o. s ] are \H\-valued and belong 
to Domo\ Moreover, (compare with Provosition \3.1\) 



E f F (B T )dB T f G (B a )dB a = a>H f f \\t — a\ 2H 2 M (r, ct)+7 (t, a) P (r, cr)l drdtr 
Jo Jo Jo Jo L J 

(3.10) 

where M(r,a) — EF (B T ) G (B a ) and 7(7", <r) (defined in (|3.4|) ) are i/ie same as m 
Proposition \S '. 11 but now (denoting by fu,V the density of a random vector (U,V)) 



P(r )( r)=f8G'(/B T , B ,):= fB r ,B„(x,y)F' (dx) G' (dy) . (3.11) 



Example 3.7. Let ,F(a;) = G(a;) = sgn(a;). Since F' = G' = 2<5 , the terms in (f37TO|) are 
given by (|3.4j) and 

F(r,a) = 4/ Br , Bg (0,0)= - - * _— , (3.12) 

irT H a H y/1 - p 2 (r, ct) 



2 r i 

M(t,<t) = -arccos V 1 ~ P 2 ( r , °") ■ ( 3 - 13 ) 



7T 



where 



R H {t,<j) /qia\ 
P\ T ^) = H— H- • ( 3 - 14 ) 



(The standard formula for the first quadrant Gaussian measure was used to ob- 
tain (|3.13[1 .) The resulting double integral cannot apparently be computed explicitly, 
but it is easy to see by substitution that 

' E I sgnB T dB T / sgnB^ dB a ^ --^ =a H \t- s\ 



dtds V ,/ n " b ' ~ ' / n b ° ~ ° / ^ Sids 
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Thus, unlike the case where H = 1 /2, X t = J^sgnB T dB T is not itself a fractional 
Brownian motion with parameter H or (as can easily be verified in the same way) 
with any other parameter H. As mentioned in the Introduction this conclusion had 
been reached in [2] via chaos expansion techniques. 



Proof of Theorem \3.6[ Assume without loss of generality that both F and G are in- 
creasing functions. Define u l n = F n (B.) l[o,t], where {F n }^ =1 is a sequence whose 
existence is assured by Lemma l3~5l To see that (ti^) n= i converges to u l in 1? {P\%), 
note that — ¥ u pointwise as n —¥ oo, and that, by (|3.8j) . 



« (r) - u* (r)) («* (a) - M 4 (a)) |a - r\ 2H ~ 2 <M 2 (1 + |fl T |) (1 + \B a \) \a - r\ 2h 1 

which is integrable in [0,T] 2 . Thus, by the dominated convergence theorem (Lebesgue 
measure), lim^^oo — = 0. Since (again by 



Vn-u^t <a H M 2 // (l + \B T \)(l + \B a \)\T-a\ 2H - 2 dr da € L 1 (P) , 
JJ[a : T} 2 

it follows from the dominated convergence theorem that lim n _>. 0o E ||tt^ — = 0. 

The next step is to show that (5u t n )'^L 1 is a Cauchy sequence in L 2 (P). For this 
purpose note that by formula (|3.3| , 

E |<5< - <5z4| 2 =£ |K - u^||^ + 

+ oh // 7 (r,<7) £ (B T ) - (S T )) « - F' m [B a ))\ drda. 

Applying Fubini's theorem (easily justifiable since the derivatives Fj, are bounded) and 
denoting by / {x, y; t, a) the density function of [B T , B a ), it remains to show that 



71,771— >CX> 



I im __ / / K (*) ~ (*)) « (tf) - (V)) i*'* (*> 2/) dxdy = 0, (3.15) 

where 



Note that 



and that 



L Lt (x, y) = a H 1 1 f (x, y; r, a) 7 (r, a) dr da. 
'[o,t] 2 



f(x,V,T,a)<f(0,0;T,a) = 



2-kt h <j h y/1 - p 2 (t, a) 

7 (r,<r) 



— dr do- < 00 (3.16) 

A) Jo r H cr ff y/1 - p 2 (r, cr) 

(the latter can be verified, using the homogeneity of the integrand, by carrying out 
an asymptotic analysis near the singularity lines). Thus the non- negative function 
L*'* (x, y) is integrable, continuous on R 2 and vanishes at infinity (this follows by 
dominated convergence from the fact that the integrand has these same properties, 
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as a function of x and y, for almost every r, a). To show (|3.15[) . L t ' t (x,y) will be 
approximated with respect to the following norm on Co^HL 1 (R 2 ): 

h\\ = IMIoo + SU P / Ifffos/)! dy + sup / |p(a;,y)| dx+ Ugl^. 

i£t Ji yeRJE. 

The set of functions which can be written as 
I 

(a;) e Co (K) n L 1 (M) , i=l,...,l, (3.17) 



(=i 



is dense (with respect to the above norm) in Co (K 2 ) fl (R 2 ) (using bump functions 
for the integer lattice squares). If, for e > 0, g e denotes such a function satisfying 
l|£ M - Sell < e , then by (f3~§)) the limit (|3~T5j) holds for gr e instead of L*»*. Moreover, 



|K (*) - K (x)| IK (y) - K (y)| IL 4 '* (z,y) - g e (x,y)\ dxdy 



< 



+ 



Fl (x) - Fl (x) Fl (y) - Fl (y) \L*>* (x, y) - g e (x, y)\ dxdy 
F n (*) - F 2 ' (x) | \Fl (y) - Fl (y) | \L^ (x, y) - g e (x, y) \ dx dy 
Fl (x) - Fl (x) Fl (y) - Fl (y) \L^ (x, y) - g t (x, y) \ dx dy 



Fl (x) - F 2 m (x) Fl (y) - Fl (y) \L^ (x, y) - g e (x, y) \dxdy< 4M 2 e 

(3.18) 

This proves (|3. 15[) . which means that (fo^J^Lj is a Cauchy sequence in L 2 (P). Since 
6 is a closed operator, it follows that u* € Dom (S) and that 

5u l = lira Sul in L 2 (P) . 

n— >oo 



Now let G n be a sequence assured by Lemma 13.51 for the function G. All the above 
holds for = G n (B.) lr 0)S i as well, and 



E5u t 5v s = lim ESu*Sv' ^ 

71— >O0 



an lim 

n— >oo 



M (r, cr) \r — a\ 2H 2 dr da 

F^x) G' n (y)L t ' s (x,y) dxdy, 



where 

L*''(x,y) 
It remains to check that 



nl(T,(j) f (x,y,T,a) drdcr. 



lim 

n— >oo 



/5 2 



K (re) G' n (y) (x, y) dxdy= / L Ls (x, y) F' ( dx) G' ( dy) . (3.19) 
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With the same argument as for L*' 4 above, one can conclude that L t,s £ Cq (IR 2 ) H 
L 1 (R 2 ) . Let g e be a function of the form (|3~T71) such that ||i 4 ' s - # e || < e. Then (|3"T31) 
holds for g e by Lemma 13.51 Denoting 



max / 


dF' 




dG' 




dA 




dA 


) 

OO 



? ^dicrete ' 



1 > G dicrete ( R ) 



it follows that 

K (x) G' n (y) L t,s (x, y) dxdy- ( L t,s (x, y) F' ( dx) G' ( dy) 

Jm 2 

■ II \F^(x)G' n (y)\\L^(x,y)-g e (x,y)dxdy\ + 
F' n (x) G' n (y) g e (x, y) dx dy - F> ® G' (g e ) 



+ \F'®G' (L*' s ) -F' ®G' (g e )\ 



<4M z e- 



F^ (x) G' n (y) g t (x, y) dxdy - F 1 ® G' (g e ) 



AK 2 e n ^f A(M 2 + K 2 )e 



where the first bound is justified exactly as p,18j) was. This proves (13.19|) and con- 
cludes the proof. □ 



4 Concluding Remarks 

We first provide an intuitive explanation of why the restriction H > x fi was neces- 
sary. As H decreases, both the integrand and the integrator in J Q * F{B T ) dB T become 
rougher, whereas the stochastic integral's existence requires a minimal amount of "cu- 
mulative" regularity among both terms (this can be best seen in Zahle's definition [S] 
of the integral). The threshold turns out to occur at H = 1 /2 (Brownian motion). 
When H < 1/2 the process B t itself does not belong to the domain of the divergence 
operator S, nor is it even an H-valued random variable. In [1) Cheridito and Nu- 
alart have extended the divergence operator to a larger domain, which includes Bt 
itself and some functions of it. However, many of the formulae and theorems don't 
carry through, and those which do take a different form. The case H < i is therefore 
significantly different. 

Secondly, and as H decreases to Y2, one expects the cross-covariance of the fractional 
stochastic integrals to converge to that of the Brownian stochastic integrals, and we 
now proceed to check, without too many details, that this indeed turns out to be the 
case. Namely, the right hand side of (|3.10[) converges to J^ M Ey 2 F(B T )G(B T ) dr (here 
and henceforth in this section, we add the underlying Hurst parameter as a subscript 
wherever relevant). For simplicity, consider first the case of Example 13.71 

i) lmiy_,.i/ 2 + an \t— <j\ 2H ~ 2 — 5(t—u) in the sense of distributions, and it is then quite 
easy to conclude that the first term in the right hand side of (|3.10|) converges to 
/; At a /2 (F(M/ r )G(^ r ))dr. 
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ii) lim aH7H(T, a) = a.e. 
ff->7 2+ 

w) lim Ph(t,(j) = 



H->i/z+ tx^Jt A cr(r V fi-r A cr) 

Items ii) and m) imply that the second integrand in (|3.10[) converges to a.e. It 
remains to find an integrable upper bound for / , uniformly for all 

T H a H y/l-~p 2 H (T,a) 

Hg [!/2, 3 /4], in order to conclude, by dominated convergence, that the integral of the 
second term converges to 0. 

iv) Note first that an"fH(T, cr) is globally bounded. Secondly, we pass to the polar 
coordinates (r, 9) for (t, <t), and split the domain t>ct>0 into the two regions 
{9G (0, ^)} and {9g , fj}. It was then possible to obtain the estimate 



1 < cfr^Vr" 



if fie (0,5) 



which, after multiplying by the Jacobian r, is in L 1 ((0, j) x (0, i?)) for any i?. 



The general case, when F and G are arbitrary elements of G, follows from this par- 
ticular example after realizing that there exist two positive constants C\ and C2 such 
that 

|F'(dx)®G'(dt/)(/B T , B J| < CxH/b^JL + O, 

= g i + ^ 

AnT H cr H y/1 - p 2 H (r, cr) 

A Appendix 

Proof of Lemma [27B . Let us first prove the "in that case" part of the proposition. First 
assume u — Fk where F is a smooth random variable and k G H. Then by definition, 
Du = DF Cg) fc and thus for each h £ H we have that (u, /i) = (fc, /i) F is a smooth 
random variable and 

L>m (/i) = (fc, h) DF = D ((k, h)F)=D ((Fk, h)) = D (u, h) 

as required. It is therefore clearly true for sums of such elements, which form a dense 
subset of D 1,2 (H). It remains to prove that if u n — > u and Du n — !• Du in I? (il, H) 
and L 2 (fi, H <g> H) respectively, such that the proposition is true for all {«„}™ =1 , then 
it is true for u. 

First note that (u n ,h) — > (u,h). To show that (u, h) € ED 1 ' 2 it will therefore suffice 
to show that {D (u n , /i)}?Lj is a convergent sequence in L 2 (f2, H) (since D is a closed 
operator, it will show that (u, h) £ D 1,2 and D(u,h) — lim^oo D (u n , h) ) . Since 
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D(u n ,h) = Du n {h), it will suffice to show that Du n (h) — ► Du(h) in L 2 (Q,H). 
Indeed, 

E \\Du n (h) - Du {h)\\ 2 =E\\D K - u) {h)\\ 2 < E\\D(u n - u)\\ 2 \\h\\ 2 



\ 2 E\\Du n - Duf 







as required. 

Now let us prove the equivalence. 

For the first implication, assume u G D 1 ' 2 (K) and let {e n }^ =1 be any orthonormal 
basis in K. Then (since the last part of the theorem is already proved - using only the 



assumption u € 



and therefore 



ill, 2 



\\Du\\ 2 = J2\\Du(e n )\\ 2 = J2\\D(u,e n ) 



n=l 



n=l 



E (J2\\D(u,e n )\\ 2 ) =E\\Du\\ 2 <^ 



as Du G L 2 (P, K (g) H). The Fubini theorem now finishes the job. 

For the second implication, assume the sum in the proposition is finite. Let 



fc=i 



for an orthonormal basis {e^}^^ of K . Then 

E\\u-u n \\ 2 , 

which converges to zero as the tail of the convergent series 



E^ 



\u\\ 



L 2 (P,H) - 



fe=l 



To show that u € D 1 ' 2 (K) it will now suffice to prove that {Du n }^ =1 converges in 
L 2 (P,K ® H). For this, let us show that it is a Cauchy sequence. Assume k > n. 
Then 



\\(Du k - Du n ) (k)\\ 2 = \\D{u k -u n ,k)\\ 2 

k 

E (e m ,k)D(u, e m ) 



D E ( U ' e m)( e m,k) 



m— n+1 



m— n+1 



<( E \(e m ,k)\\\D(u,e m )\\\ 

\m— n+1 / 



< E i( e - fc >i 2 E h d < 



11. c Tl 



\m—n-\-l 



\m— n+1 
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whence 

E\\(Du k -Du n )(k)\\ 2 <( J2 \( e m,k)A ( J2 E\\D(u,e m }\\ 2 

\m— n+l / \m— n+1 

with both sums being finite. It is now clear that the sequence in question is indeed a 
Cauchy sequence, whence we are done. □ 
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